Knudsen Diffusion in Silicon Nanochannels 
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Measurements on helium and argon gas flow through an array of parallel, linear channels of 12 nm 
diameter and 200 fim length in a single crystalline silicon membrane reveal a Knudsen diffusion type 
transport from 10 2 to 10 7 in Knudsen number Kn. The classic scaling prediction for the transport 
diffusion coefficient on temperature and mass of diffusing species, D hc oc \/r i s confirmed over a 
T-range from 40 K to 300 K for He and for the ratio of Due/DAi oc ^JnriAr /mn c . Deviations of the 
channels from a cylindrical form, resolved with electron microscopy down to subnanometer scales, 
quantitatively account for a reduced diffusivity as compared to Knudsen diffusion in ideal tubular 
channels. The membrane permeation experiments are described over 10 orders of magnitude in Kn, 
encompassing the transition flow regime, by the unified flow model of Beskok and Karniadakis. 

PACS numbers: 47.45.n, 47.61.-k, 47.56. +r 



Knudsen diffusion (KD) refers to a gas transport 
regime where the mean free path between particle- 
particle collisions A is significantly larger than at least 
one characteristic spatial dimension d of the system con- 
sidered By virtue of the negligible mutual particle 
collisions transport in such systems takes place via a se- 
ries of free flights and statistical flight direction changes 
after collisions with the confming walls. Because of the 
dependency of A on pressure p and temperature T given 
by the kinetic theory of gases, Å oc T/p, this transport 
regime is only observable in macroscopic systems at very 
low p or elevated T. By contrast, for transport in spa- 
tially nanoconfined systems with d = O(10 nm), the 
Knudsen number Kn — X/d, which quantifies the gas 
rarefaction, is larger than 1 even for ambient pressures 
and temperatures, e.g. A(He) = 118 nm at p — 1 bar 
and T — 297 K. Thus, for many processes involving 
gas transport in restricted geometries, e.g. gas catalysis 
and storage or equilibration phenomena via gas flow in 
meso- and nanoporesJH KD plays a crucial role. 

Albeit the phenomenon of KD has been known for al- 
most a hundred years, details on its microscopic mecha- 
nisms are still at debate, in particular how wall rough- 
ness down to atomic scales and particle-wall interaction 
parameters affect the statistics of the diffusion process 
and thus the value of the KD transport coefficient D [B| . 
In the past, a comparison between theory and experi- 
ment has often been hampered by complex pore network 
topologies, tortuous transport paths and poorly charac- 
terized pore wall structures in available nanoporous ma- 
trices. Nowadays, membranes with better dermed geome- 
tries, such as carbon nanotube bundles [^Ø], porous alu- 
mina 8] , and tailored pores in mesoporous silicon 0, Q , 
allow us to address such fundamental questions in more 
detail. A better understanding of, and more in-depth in- 
formation on, transport in such systems is not only of 
academic interest, but also of importance for the archi- 
tecture and functionality in the emerging held of nano- 




electromechanical and nanofluidic systems 11 1 . 



FIG. 1: (color online) He pressure relaxations in RI and R2 
at T — 297 K for starting pressures p s — 0.1 mbar (stars) 
and 100 mbar (triangles) without (a) and with SiNC mem- 
brane (b). (c) Pressure relaxation time r as determined from 
measurements without (squares) and with SiNC membrane 
(circles) versus Knudsen number in the capillaries Kn c (bot- 
tom axis) and in the SiNCs Kn (top axis). The lines in (a), 
(b) and (c) represent calculated p(t) and r(Kri) values. Insets 
in (c): TEM cross section views of two SiNCs in comparison 
with circular perimeters of 12 nm diameter (dashed lines). 



In this Letter, we report the first rigorous experimen- 
tal study on transport of simple rare gases, i.e. He and 
Ar through an array of parallel aligned, linear channels 
with ~ 12 nm diameter and 200 /im length in single crys- 
talline Si over a wide temperature (40 K < T < 300 K) 
and Kn-r&nge (10 2 < Kn < 10 7 ). We explore the clas- 
sical scaling predictions given by Knudsen for D on T, 
mass of diffusing species, m, and explore its dependency 



on Kn. Attention shall also be paid to the morphology 
of the channel walls, resolved with transmission electron 
microscopy (TEM), and how it affects the diffusion. 

Linear, non-interconnected channels oriented along the 
(100) Si crystallographic direction (SiNCs) are electro- 
chemically etched into the surface of a Si (100) wafer 
[il]. After the channels reach the desired length L of 
200 ± 5 fim, the anodization current is increased by a 
factor of 10 with the result that the SiNC array is de- 
tached from the bulk wafer [13| . The crystalline Si walls 
call for irregular channel perimeters as one can see in 
the TEM pictures of Fig. QJc), recorded from a mem- 
brane part Ar ion millcd to ~ 5 nm thickness. In first 
approximation, however, the SiNC's cross-section can be 
described as circular with a diameter of ~ 12 nm, in ac- 
cordance with the analysis of an Ar sorption experiment 
at T = 86 K, which yields additionally a SiNCs density 
of $ = (1.5 ±0.1) • 10 11 cm" 2 . 

Our experimental setup consists of a copper cell with 
an inlet and outlet opening 14|. Inlet and outlet are 



connected via stainless steel capillaries of radius r c = 
0.7 mm and length l c — 70 cm with two gas reservoirs, 
RI and R2, of an all-metal gas handling. Additionally 
two pneumatic valves VI and V2 are used to open and 
close the connections between the sample cell and RI and 
the sample cell and R2. Four thermostatted capacitive 
pressure gauges allow us to measure the gas pressures 
in RI and R2, p\ and p2, resp., over a wide pressure 
range (5T0 -3 mbar < p <1 bar) with an accuracy of 
10~ 3 mbar. The cell is mounted in a closed-cycle He 
cryostat to control the temperature between 40 K and 
300 K with an accuracy of 1 mK. 

Our goal of studying gas transport dynamics over a 
large Kn range necessitates a thorough understanding of 
the intrinsic flow characteristics of our apparatus. There- 
fore, we start with measurements on He flow through 
the macroscopic capillaries and the empty sample cell at 
room temperature T — 297 K. We record the equilibra- 
tion of pi (t) and P2 (t) towards a pressure p eq as a function 
of time t after initial conditions, p\ — p s > pi — mbar 
at t = s. In Fig.[T](a) px(t)/p s andp 2 (i)/p s are plotted 
for starting pressures of p s — 0.1 and 100 mbar. We ob- 
serve pressure equilibrations towards p cq — 0.4 ■ p s . The 
value 0.4 is dictated by the volume ratio of RI to R2. 
From the p(i)-curves we derive characteristic relaxation 
times r according to the recipe p\ (t = t) — pi (t — r) = 
l/10-p s . It is understood that p is monotonically decreas- 
ing downstream from RI to R2. In order to quantify the 
gas rarefaction we resort, therefore, to a calculation of a 
mean Knudsen number Kn c — X(p)/r c in the capillaries 
assuming a mean pressure, p — p cq . This simplification 
is justified by the analysis provided below and by theo- 
rctical studies which indicate differences of less than 1% 
between flow rates calculated with an exact and an av- 
eraged treatment for Kn lfj, 15j . In Fig. [1] (c) we plot 



ation from 0.01 to 100 mbar. For Kn c < 0.1, r increases 
with increasing Kn c . In an intermediate Kn c range, 
0.1 < Kn c < 0.6, we observe a cross-over regime towards 
a saturation plateau with r ~ 140 s that extends to the 
largest Kn c studied. These changes in the p relaxation 
and hence flow dynamics are reminiscent of the three 
distinct transport regimes known to occur for gases as a 



function of thcir rarefaction 



lffl: For Kn c < 0.1 the 



number of interparticle collisions still predominates over 
the number of particle-wall collisions. Hagen-Poiscuillc's 
law is valid and predicts a decreasing flow rate and hence 
increasing r due to the 1 /ifn-scaling of the particle num- 
ber density in gas flows. For Kn c > 1 we enter the pure 
KD regime|l], wherc theory predicts a He KD transport 
coefficient Df c dependent, however a Kn c independent 
particle flow rate, 
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which is responsible for the plateau in r for the larger 
Kn c investigated. In Eq. [T] p- u p Q refers to the inlet 
and outlet pressure of the capillary considered, resp., 
and &b to the Boltzmann factor. In the intermediate 
Kn-r&nge, the interparticle collisions occur as often as 
particle-wall collisions which gives rise to the cross-over 
behavior found for r. 

We elucidate this behavior in more detail by divid- 
ing the flow path within our apparatus into two flow 
segments (up- and downstream capillary) and calculate 
the particle number changes along the flow path and 
the resulting pi(t), P2{t) with a 1 ms resolution using 
the unifiedflow model of Beskok and Karniadakis (BK- 
model) |T(J E3, El and a local K n-number for each flow 
segment, Kn\ = Kn({p\ +p Q )/2): 
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measured r values versus Kn c corresponding to a p s vari- 



Eq. [2] comprises a Hagen- Poiseuille term, a term which 
treats the transition of the transport coefficient from 
continuum-likc, i. c. the bare dynamic viscosity n(T), 
to the KD transport coefficient, D^ c with a = 
2.716/7rtan _1 (aiifn' 5 ), and a generalized velocity slip 
term, which is second-order accurate in Kn in the slip 
and early transition flow regimes (Kn < 0.5). The model 
captures for Kn — > the no-slip Hagen-Poiseuille limit, 
whereas it transforms to Eq. Q] for Kn ^> 1. As veri- 
fied by comparison of the BK-model with Direct Simula- 
tion Monte Carlo and solutions of the Boltzmann equa- 
tion, a choice of fe = —1 for the slip parameter results 
in the correct velocity profile and flowrate, as well as a 
proper pressure and shear stress distribution in a wide 
-Kn-range, including the transition flow regime. After 
optimizing the free parameters in Eq. Ø a± and /?, we ar- 
rive at the p(t)- and r(-ftTn c )-curves depicted in Fig.[T](a) 
and (c), resp. The good agreement of the BK-model pre- 



dictions with our measurements is evident and the ex- 
tracted parameters a.\ = 10 ± 0.2 and (3 = 0.5 ± 0.05 
(Z)He _ o.412 ± 0.06 m 2 /s) agree with BK-modelling of 



He gas flow as a fvmction of rarefaction 
ogous analysis for Ar gas flow yields ol\ 
$ = 0.1 ± 0.02 {D^ = 0.134 ± 0.02 m 2 /s). 
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FIG. 2: (color online) Pressure relaxation time r for He (cir- 
cles) and Ar (triangles) measured with built-in SiNC mem- 
brane at T = 297 K versus Knudsen number in the SiNCs Kn 
(bottom axis) and in the capillaries Kn c (top axis). The lines 
represent calculated r(7fn)-curves. Inset: TEM cross section 
of a SiNC in comparison with a circular channel perimeter 
of 12 nm diameter (dashed line). The SiNC's perimeter is 
highlighted by a line. 

Having achieved a detailed understanding of the intrin- 
sic flow characteristics of our apparatus, we can now turn 
to measurements with the SiNC membrane. The mem- 



brane is epoxy-sealed in a copper ring 14] and special at- 
tention is paid to a careful determination of the accessible 
membrane area, A — 0.79 ± 0.016 cm 2 in order to allow 
for a reliable determination of the number of SiNCs in- 
serted into the flow path, N = $-A = (11.8±1)-10 10 . As 
expected the pressure equilibration is significantly slowed 
down after installing the membrane - compare panel (a) 
and (b) in Fig. [1] Choosing selected p s within the range 
0.005 mbar to 100 mbar, which corresponds to a variation 
of Kn in the SiNCs from 10 2 to 10 7 , we find an increase 
in t of ~ 170 s, when compared to the empty cell mea- 
surements - see Fig. [T] (c). From the bare offset in r for 
Kn c > 1, we could calculate Dn c in the SiNCs. We are, 
however, interested in the Due behavior in a wide Kn- 
range, therefore we modify our flow model by inserting 
a segment with KD transport mechanism charactcristic 
of N tubular channels with d = 12 nm and L — 200 /im 
in between the capillary flow segments. Adjusting the 
single free parameter in our simulation, the value of the 
He diffusion coefficient Z?h c in a single SiNC, we arrive 
at the p(t)- and r(i4Tn)-curves presented as solid lines in 
Fig. [T] (b) and (c), resp. The agreement with our mea- 
surement is excellent and the model yields a Kn inde- 
pendent D Ho = 3.76 ± 0.8 mm 2 /s. 



It is worthwhile to compare this value with Knudsen's 
prediction for D. In his seminal paper he derives an ex- 
pression for D with two contributions, a factor G, charac- 
teristic of the KD effectivity of the channel's shape, and 
a factor solely determined by the mean thermal velocity 
of the particles v given by the kinetic theory of gases, 



D=\Gv=\G 
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Interestingly, by an analysis of the number of particles 
crossing a given section of a channel in unit time after 
complctcly diffuse reflection from an arbitrary element of 
wall surface and while assuming an equal collision acces- 
sibility of all surface elements, Knudsen derived an ana- 
lytical expression for G for a channel of arbitrary shape. 
Knudsen's second assumption is, however, only strictly 
valid for circular channel cross-sections, as pointed out 
by v. Smoluchowski [l[ and elaborated for fractal pore 
wall morphologies by Coppens and Froment Given 
the roughly circular SiNCs' cross section shapes, we nev- 
crthcless resort to Knudsen's formula, here quoted nor- 
malized to the KD form factor for a perfect cylinder: 
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The integral in Eq. [4] depends on the ratio of perimeter 
length o(l) and cross sectional area A(l) along the chan- 
nel's long axis dircction l, only. This ratio is optimizcd 
by a circle, accordingly the most effective KD channel 
shape is a cylinder, provided one calls for a fixed A along 
the channel. Eq. |4] yields due to our normalization just 
G = d and we would expect a Z?Hc of 5 mm 2 /s for He KD 
in a SiNC, if it were a cylindrical channel of d = 12 nm. 
A value which is ~ 30% larger than our measured one. 
If one recalls the TEM pictures, which clearly indicate 
non-circular cross-sections this finding is not surprising. 
In faet these pictures deliver precisely the information 
needed for an estimation of the influence of the SiNC's 
irregularities on the KD dynamics. We determine the ra- 
tio o/A of 20 SiNC cross sections, and therefrom values 
of G. The values of the SiNC in Fig. Q] (top), (bottom) 
and Fig. [2] corresponds to Gs of 9.1 nm, 10.3 nm, and 
10.6 nm, resp. Tacitly assuming that the irregularities 
exhibited on the perimeter are of similar type as along 
the channel's long axis we take an ensemble average over 
the 20 G values and arrive at a mean G of 9.9 nm, which 
yields a Z?Hc of 4.1 mm 2 /s, a value which agrees within 
the error margins with our measured one. 

To further explore the KD transport dynamics in the 
SiNCs we now focus on the m and T dependency of 
D. In Fig. [2] the r(Kn) curve recorded for Ar and He 
at T = 297 K are presented. Both exhibit a similar 
form, the one of Ar is, however, shifted up markedly 
towards larger t. Our computer model can quantita- 
tively account for this slowed-down dynamics by a fac- 
tor 2.97 ± 0.25 decrease in D as compared to the He 



measurements, which confirms the prediction of Eq. [3l 
Dn e /DAr = ^/niAr/wHc = 3.16. For the exploration of 
the T behavior we choose again He as working fluid due to 
its negligible physisorption tendency, even at low T. We 
perform measurements at selected Ts from 297 K down to 
40 K, shown in Fig.O We again perform computer calcu- 
lations assuming KD in the SiNC array superimposed to 
the transport in the supply capillarics, prcscntcd as lines 
in Fig.[3j and optimize the single frec parameter D(T) in 
order to match the observed r(Kn, T) behavior. Despite 
unresolvable deviations at larger Kn c and decreasing T, 
which we attribute to thermal creep effects, characteristic 
of T gradients along the flow path [1(3] , we find, in agree- 
ment with the experiment, increasingly faster dynamics 
with decreasing T, Note, this counter-intuitive Ånding 
for a diffusion process results from the l/T scaling in 
Eq. Q] which reflects the T-dependency of the particle 
numbcr density in gas flows. By contrast, the Dn c (T), 
determined by our simulations, scales in excellent agree- 
ment with y/T, see inset in Fig. [3j confirming Knudsen's 
classic result down to the lowest T investigated. 
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FIG. 3: (color online) Pressure relaxation time r for He gas 
flow through SiNC membrane (symbols) and simulated r(Kn) 
(lines) versus Knudsen number in the SiNC membrane Kn 
(bottom axis) and in the supply capillaries Kn c (top axis) at 
selected temperatures T. Inset: He diffusion coefflcient Due 
in SiNCs in comparison with the y/T prediction of Eq.[3](lme) 
plotted versus y/T. 

We find no hints of anomalous fast KD here, as was re- 
cently reported for a bundle of linear carbon nanotubes 
and explained by an highly increased fraction of spec- 
ularly reflected particles upon wall collisions [l[ . The al- 
tered collision statistics was attributed to the crystalline 
structure and the atomical smoothness of the nanotube 
walls. The SiNC walls are also crystalline, however, not 
atomical flat, as can be seen from our TEM analysis. 
Along with the formation of a native oxide layer, typical 
of Si surfaces, which renders the near surface structure 
amorphous, silica like, this, presumably, accounts for the 
normal KD observed here contrasting the one found for 
the graphitic walls of carbon nanotubes. 

We presented here the first detailed study of gas trans- 



port in linear SiNCs. Our conclusions are drawn from a 
correct treatment of gas flow over 10 orders of magnitude 
in gas rarefaction, which is, to the best of our knowledge, 
the largest Kn range ever explored experimentally. The 
characteristic properties of KD, an independency of D 
on Kn, its scaling predictions on to, T, and on details of 
the channel's structure, resolved with sub-nm resolution, 
are clearly exhibited by our measurements. 
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